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ABSTRACT 


Techniques  for  determination  of  "dynamic"  error  parameters,  or  those 
obtained  from  real  trajectories,  are  studied  in  detail,  with  particular 
attention  being  given  to  dependence  of  estimates  on  trajectory  ar.d  to 
shortcomings  of  est  mation  procedures  which  ignore  such  dependence. 

A  technique  is  presented  to  determine  dynamic  estimates  of  error  para¬ 
meters  for  the  rawiniionde  and  Is  verified  by  data.  The  data  indicate 
that  error  paramete-s  so  determined  are  smaller  than  given  by  most 
earlier  studies  and  tend  to  conform  to  error  ciaims  made  for  the 
AN/GMD-UA)  Raw  inset. 
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INTRODUCTION 

If  a  moving  object  is  under  observation  by  a  given  tracking  system,  the 
observations  so  made  generally  take  the  form  of  a  sequence  of  vectors, 
each  of  which  is  related  to  a  specific  time,  and  whose  components  give 
estimates  of  tne  object’s  position  in  some  convenient  coordinate  system. 
For  a  radar,  this  sequence  is  of  the  form  {(Aj,  £j,  e.j)|i=l,n},  where 
each  vector  ir.  the  sequence  represents  range,  azimuth,  and  elevaticr,. 
respectively,  at  times  {tj|i=l,n}  in  the  obvious  fashion.  In  contrast, 
the  rawinsonde  system  fqr  balloon  tracking  gives  rise  to  the  family 
{(a'j,  e’j,  £}..  p j,  fej)|i=I,n},  where  the  components  of  a  vector  re¬ 
present  azimu+n  and  elevation  from  the  AN/GMD-(  )  windset,  and  temper¬ 
ature,  pressure  and  humidity  from  the  radiosonde  balloon.  Here  agcin 
all  vectors  are  indexed  by  a  time  value.  From  each  vector  (a'j,  e'j, 

£' j,  pj,  h\) ,  a  new  vector  (a'j,  Zj)  is  derived  through  use  of 

the  hydrostatic  equation.  The  family  {(a’j,  e.' j ,  z  j )  1 1  =  i ,  n>  so  formed 
represents  position  coordinates  in  a  curvilinear  system.  It  is  this 
system  which  will  be  considered  in  the  following  study. 

In  the  statistical  model  to  be  employed,  a  sequence  of  vectors  arising 
in  the  fashion  aiscussid  previously  will  be  considered  as  a  single 
realization  of  an  appropriate  subprocess  of  a  multivariate  stochastic 
process.  The  index  set  for  this  process  will  be  a  set  T  containing 
(t j | i= I , n}.  For  each  t  in  T,  the  corresponding  multivariate  random 
variable  will  have  as  its  mean  a  vector  which  represents  rhe  position 
of  the  object  in  an  appropriate  coordinate  system  and  wiL  have  as 
its  variance-covariance  matrix  a  matrix  Z  independent  of  t.  The  pro¬ 
blems  to  be  addressed  in  this  study  are  concerned  with  obtaining  es¬ 
timates  of  the  matrix  Z  from  a  single  realization  of  some  multivariate 
subprocess.  Such  estimates  will  be  designated  as  "dynamic",  as  they 
represent  the  behavior  of  the  system  in  motion,  as  contrasted  to 
slatic  tests  which  obtain  error  estimates  for  a  fixed  configuration. 

The  marginal  processes  are  in  general  nonstationary,  even  in  the  weak 
sense  CG  due  to  the  variation  of  mean  values  with  changing  t,  so 
that  the  problems  involved  may  be  considerable. 

Chief  among  the  problems  encountered  is  the  dependence  of  the  estimates 
on  trajectory.  SuopGse  it  is  desired  to  obtain  dynamic  error  parameters 
for  a  radar.  Superficially,  one  might  adopt  the  following  procedure. 

The  radar  to  be  tested  and  a  second  radar  of  known  error  behavior  are 
placed  in  close  proximity  and  allowed  to  track  an  object  simultaneously 
for  a  time  interval  T.  The  radars  may  be  expected  to  g.ve  rise  to 
vector  sequences  {(/tj,  o.j ,  ej)|i=l,n}  and  {(A'j,  a'j,  e'j)|i=l,n}, 
respectively,  where  the  unprimed  sequence  will  represent  the  radar  of 
unknown  error  characteristics.  As  position  and  velocity  are  generally 
studied  in  a  cartesian  coordinate  system,  transform  each  to  a  single 
fixed  reference  system  (X,  Y,  Z),  obtaining  sequences  {(Xj,  y j,  Zj)|i=l,n} 


( 


and  {(x* | ,  j ,  z * i ,-> 1 1=1, n}.  If  calibration  procedures  are  correct,  up 
to  errors  introduced  by  the  system  itself,  these  values  should  be  . 
identical.  Therefore,  members  of  the  vector  sequence  { (Xj— xT j ,  y\-y'\, 
Zj-z? j ) | i=| ,n}  should  have  near-zero  components,  and  the  variances  and 
covariances  required  to  determine  the  error  behavior  of  the  radar  ,with< 
reference  to  (x,  y,  z>  position  estimates  may  be  calculated.  With 
these  values,  the  behavior  of  derivative  estimates  are  determined.  This 
procedure  has  been  in  common  use  in  the  past  to  estimate  error  behavior. 

To  understand  the  difficul+ies  introduced  by  such  a  procedure,  additional 
background  will  be  required  and  will  be  presented  in  the  following  section 
This  presentation  will  be  made  for  two  reasons.  First,  it  is  worthwhile 
from  the  standpoint  of  understanding  error  behavior  introduced  by  mathe¬ 
matical'  transformations  in  a  radar,  and  second  it  is  necessary  for  the 
understanding  of  a  technique  to  be  introduced  for  determining  dynamic 
error  parameters  for  the  raw  insonde.  ' 


ERROR  BEHAVIOR  IN  THE  TRANSFORMATION  FROM 
SPHERICAL  TO  RECTANGULAR  CARTESIAN  COORDINATES 

Let  an  object  be  moving  along  a  trajectory  {(p(t),  a(t),  e(f))|teT)  j 
relative  to  a  given  spherical  coordinate  system.  Suppose  a  radar'  is 
located  at  the  origin  of  this  coordinate  sys+em  and  oriented  in  such  a 
fashion  that  a  point  (p,  a,  e)  in  the  system  corresponds  to  range, 
azimuth,  and  elevation  as  seen  by  the  radar.  Let  a  sequence 
{(/Lj,  £j)|i»|,n}  of  observations  be  taken  in  the  time  interval  T, 

along  the  given  trajectory,  these  observations  occurring  at  times 
{t j  j i= I , n> .  It  is  desired  to  determine  the  trajectory  of  the  noving 
object  with  respect  to  a  rectangular  coordinate  system  related  to  the 
first  by  the  transformation  .  . 

/ 

x  =  pcosacose  : 

/y  =  psinacose  ( | )' 

z  =  psine 

which  wi 1 1  be  cal  led  T. 

From  (I)  the  family  ,  £j)|i  =  !,n)  gives  rise  to  a  new  family 

{(Xj,  yit  Zj)|i=i,n}  in  an  obvious  manner.  If  the  first  of  these  fami¬ 
lies  is  considered  as  a  realization  of  a  subprocess  of  a  combined 
random  process  {(R+,  A+,  E-j-)|teT},  { ( x  j ,  y  j,  Zj)|i  =  l,n}  becomes  a  reali¬ 
zation  of  a  subprocess  of  the  joint  stochastic  process  {(Xf,  Vf,  Z+)|teT} 
defined  for  each  +  in  T  by  . 


I 


1 


X,  =  R, cos A, cos E. 
t  f  ■-  f  t 

’  V^.  =  R^sinA^cosE^. 

j 

Z,  =  R.sinE.  : 

It  t  t 


Recall,  that  for  any  t  in  T,  (R-f.,  Ay,  E-j.)  has  mean  vector  (p(t),  a(t), 
c(t))  and  variance-covariance  matrix  Z.  Let  R  be  a  region  in  three-spare 
in  which  T  is  at  least  of  class  C1  (possesses  continuous  partial  deriv¬ 
atives  of  order  at>  lea'st  l)?  and  assume  {(p(t),  a:(t),  e(t))|teT}  lies 
entirely  in  R.  It  is  known  C23  that  for  a  given  precision  figure’  6, 
there  exist  a  transformation  L  and  a  constant  X  dependent  on  6  such 
that 


Pr{|lL(R+,  A+,  E+)  -  T(R+,  A,f,  E+)  1 1  H}  ^ 


tr(E) 


where  tr(£)  denotes  the  trace  of  I'. 

*  i 

From  (2),  if  tr(E)  is  sufficiently  small,  the  random  vectors  ,L(Rf,  A,,  Ef) 
and  T(R+,  Af,  Ef)  differ  only  Slightly  in  the.ir  behavior.  The  advantage 
of  this  lies  iri  the  fact  thpt  the  variance-covariance  matrix  of  URy,  Af,  Ef) 
has  the  simple  form  dT+EdT+,  where  dT+  is  the  differential  of  the  trans¬ 
formation  T  evaluated  at  (p(t),  a(t),  e(t)),  and  dT+  denotes  the  trans¬ 
pose  of  dTf  [3].  The  matrix  dT-f-  has  a  simple  form 

(cosa(t)coss(t)  -p(t)sina(t)cose(t)  -p(t)ax  a(t)sine(t)l 


dT^.  =  .  sina(t)cose(t)  p(t)cpsa(t)cose(t)  -p(t)sina(t)sine(t)  (3) 


jsine(f) 


p(t)cose(t) 


From  (3)  the  matrix  2'f  =  dTfZdT+  may  easily  be  calculated. 

As  an  example,  let  the  object  tracked  be  a  balloon  risjng  at  a  fixed  rate, 
i  a,  in  still  air.  Assume, the  ba I  loon  is  released  , a  horizontal  distance, 
b,  along  the  line  a  =  0.  The  parametric  equations  for  the  balloon  tra¬ 
jectory  are  •  '  1 

i  «  fp(t)  =  (b2  +'a  2f2)l/2 


o(t)’=  0 


e(t),  =  arcro: 


£b2  +  2*t2j 


4 


case  for  C|j  $  bc^.  Clearly  a'22(+)  remains  constant  for  all  values 
of  t  in  T  while  0 * 33  increases  from  b2a2^  at  t  =  0  approaching  c2 u 
+  b20233  as  i  becomes  large.  The  nonzero  covariance  term  a'  |3(t)=0*3i  (t> 
like  0’||(t),  achieves  an  extremal  value  only  if  d  j  j  >  b033,  in  this  case 
a  maximum,  then  decreases  without  bound  with  increasing  t  values.  In 
this  case,  when  t  is  large,  these  terms  resemble  -bato33.  If  cr||<bu33., 
a,(3(t>  =  a  *  3  j (t)  decreases  without  bound  from  t  =  0,  again  resembling 
-bara33  when  t  is  large.  It  is  not  difficult  to  see  the  role  played 
by  b  in  the  mairix  Vf.  This  emphasizes  the  dependence  of  the  matrix 
E’t  on  the  trajectory  chosen.  Clearly,  if  two  radars  at  different  loca¬ 
tions,  say  horizontal  distances  b|  and  b2,  observe  the  same  balloon, 
then  at  a  given  time  t,  the  associated  variance-covariance  matrices  of 
the  two  random  vectors  produced  would  vary  considerably.  It  would 
appear  then  that  a  realization  of  the  process  {(X-j-,  JAj.,  Z+)|teT}  is 
useless  in  the  estimation  of  error  parameters  except  in  the  grossest 
sense,  so  that  other  techniques  must  be  utilized  to  describe  the  error 
behavior  of  the  system.  Such  a  technique  will  now  be  presented. 

It  is  possible  to  approach  the  problem  from  a  direct  statistical  stand¬ 
point,  and  results  of  some  generality  may  be  obtained.  However,  these 
techniques  require  knowledge  or  approximate  knowledge  of  trajectory  and, 
for  the  case  of  the  rawinsonde,  have  the  unfortunate  characteristic  of 
being  unable  to  exclude  induced  balloon  oscillations.  For  this  case 
the  results  are  included  only  as  an  appendix. 


ESTIMATION  OF  DYNAMIC  ERROR  PARAMETERS  FOR 
THE  AN/GMD- (  )  BY  USE  OF  A  RADAR 

1 f  a  radar  such  as  the  T-9  or  FPS-16  is  available  and  in  a  proper  con¬ 
figuration  with  the  AN/GMD-(  ),  the  following  approach  to  +he  evaluation 
of  error  parameters  of  the  AN/GMD- (  )  shows  promise.  Let  a  radiosonde 
balloon  be  released  and  simultaneously  tracKed  by  the  radar  and  the 
AN/GMD-(  ).  The  output  of  the  former  of  these  devices  consists  of  a 
sequence  of  vectors  { (/tj ,  a\,  ej ) | i  =  I , n >  where  the  components  of  a  vector 
(/q,  a\,  e |)  represent  range,  azimuth,  and  elevation  estimates,  respec¬ 
tively,  for  a  time  t|  in  the  interval  C+|,  +n3  under  consideration. 

The  latter  system  has  as  its  output  the  vector  sequence  {(6|,  ^ j ) | i  =  I , n } . 
In  this  case,  the  components  of  a  vector  (b|,  ;{|)  are  estimates  of 
azimuth  and  elevation,  respectively,  at  a  time  tj  identical  to  that  time 
specified  above.  The  vector  sequences  {(/tj,  a\,  £|)|M,n}  and 
{(6j,  will  be  considered  as  realizations  of  appropriate  vector 

subprocesses  of  the  continuous  vector  processes  {(R+,  Af,  E^.)|teT}  anc 
{(8+,  F.}.)  |teT},  respectively,  where  [_+|>  +pJ  is  contained  in  T.  The 
mean  for  the  former  of  these  processes  will  be  the  vector  function 
{(p(t),  a(t),  e(t))|teT}  which  gives  the  true  trajectory  of  the  balloon 
in  the  obvious  spherical  coordinate  system  associated  *ith  the  radar. 

The  vector  function  { ( 0(  t) ,  4>{  t)  ?  |  teT}  plays  an  exactly  similar  role 
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with  the  GMD.  It  will  also  be  assumed  fhat  the  former  sequence  has 
variance-covariance  matrix  Z  independent  of  t  in  T,  and  the  latter  the 

variance-covariance  matrix  r,  likewise  independent  of  t  in  T.  Now  let 

(X,  6,  n)  be  a  vector  giving  in  component  order  range,  azimuth  and  ele¬ 
vation  or  the  radar  under  consideration  with  respect  to  the  GMD.  In 
other  words,  (X,  S,  n)  are  the  spherical  coordinates  of  the  radar  with 

respect  to  the  obvious  system  located  at  the  GMD.  (See  Figure  I.)  It 

will  be  assumed  throughout  that  the  radar  and  the  GMD  are  in  parallel 
orientation.  If  is  clear  that  if  X  is  sufficiently  large,  this  v;ou!d 
not  in  general  be  the  case,  this  due  to  the  curvature  of  the  earth. 

The  equations  relating  (3(t),  <f>(t) )  and  (p(t),  a(t),  e(t))  will  be  more 
complicated  when  the  earth’s  curvature  is  included  than  when  it  is  not; 
therefore,  from  a  standpoint  of  simplicity,  the  parallel  orientation 
model  is  desirable.  Also,  as  will  be  made  clear  later,  if  X  is  chosen 
too  large,  the  method  for  evaluating  dynamic  parameters  to  be  presented 
now  will  fail.  Then  concern  will  be  only  with  X  so  small  that  curvature 
of  the  earth  presents  a  negligible  effect. 


Under  these  circumstances,  and  under  suitable  restrictions,  the  follow¬ 
ing  equations  may  be  shown  to  hold: 


3(t)  ~  arctan 


p(t)cose(t)sina(t)  +  Xcosnsind 
p“(  t )  cose  ( t ) cosa  ( t )  +  Xcosncosfi 


d(t)  =  arctan 


p(t)sine(t)  +  Xsinn 


(4) 


where  :ji(t)  is  given  by 

i|/(t)  =  [(p(t)cose(t) ) 2  +  (Xcosn)2  +  2Xp(t)cosncose(t)cos(a(t)-6)]. 


The  equations  in  (4)  define  a  transformation  of  Euclidian  3-space  into 
Euclidian  2-space  which  will  be  denoted  by  T*.  It  follows  by  definition 
that  {T*(p(t),  a(t),  e(t))|teT}  is  Identically  the  vector  function 
{(3(t),  <J>(t))|teT}  If  quadrant  ambiguities  are  considered.  The  vector 
process  {T*(Rf,  A^,  E^.)|teT}  would  be  expected  to  exhibit  behavior  in 
its  mean  similar  to  {(8+,  Ff)jteT}.  The  complicated  nature  of  T*  pro¬ 
hibits  direct  calculation  of  the  mean  and  variance-covariance  matrices 
of  the  various  vectors  in  {T*(Rf,  Af,  Ef)|teT}  even  under  assumptions 
of  Joint  normality  of  the  process  {(R+,  Af,  E+)|teT}.  In  this  case  il 
becomes  necessary  to  resort  to  approximation.  The  magnitude  of  fhe  ele¬ 
ments  of  Z  suggests  that  satisfactory  approximations  of  both  the  mean  and 
the  variance-covariance  matrix  of  (T*(Rf,  A-j.,  E^.)|teT}  are  possible  in 
terms  of  Z  and  {(p(t),  a(t),  e(t))|teT}  f2].  In  this  light,  the  approx¬ 
imation  to  the  desired  variance-covariance  matrix  of  T*(Rf,  Af,  E^.) 
is  simply  dT*-t.Zcfir*+  where  dT*+  is  the  differential  of  the  transformation 
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FIGURE  t.  RELATED  MEASUREMENTS  OF  (0,*)  FOR  THE  A»VGMD-( 
AND  (p,a,«)  FOR  THE  AN/^PS  16  RADAR. 


AN/GMD-( 


T*  DD  a4-  (p(t),  a(t),  e(+)),  and  dT*+  denotes  the  transpose  of  this 
matrix.  The  calculation  of  dT*.j.  is  routine  but  tedious,  and  the  follow¬ 
ing  definitions  will  be  made  to  simplify  notation. 

A(t)  =  (p(t)cose(t))2  +  (Acosn)2  +  2Acosnp(t)cose(t)cos(a(t)-6) . 

For  i=l,  2,  j=l,  2,  3,  let  y*<i,j,t)  designate  the  i,  jth  entry  cf  dT*+. 
Then 


1  i  \  _  Acosncose(t)sin(a(t)-6) 

y  u,i,t;  -  ^p) 

o  Cp(t)cose(t5D2  +  p(t) Acosncose(t)cos(a(t)-S) 
y  u,z,t;  -  —  - 

_  Acosnsine(t)sin(a(t)-<S) 
y  ' 

y*l2  I  t)  - 

A/w't) 

Ap(t)cose(t)tsinnoose(t)  -  cosnslne(t)cos(a(t?--S)] 

A/^TTT 

y*(2  2  t)  -  Cp(t)slne(t)  +J  Asinn3[Ap(t)cosricose(t)sin(a(t)-6)] 

A^JTtT 

'  \ 

* 

y*(2,3,t)  =  P3(t)cos^t)  +  2Ao2(t)cos2e(t)cosncos(a(t)-6) 

Av'ijTfT 

+  M^)Cp(t)sine(t?  t  AsinnXslnncose(t)  t  cosesln£(t)cos(a(t)~6)] 

A/tnr 

Observe  fhat  the  following  relations  hold. 

I imy*( I ,2,t)  =  I imy*(2,3,t)  =  I 

p(t)->v»  p(  +  )-*o 

I  imy*(  I  ,j,t)  =  I  imy*(2',k,t)  =  0  j72,  M3. 
p(t)-H°-  p{t)-*^ 
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Hence  in  the  case  of  very  large  p(+),  the  matrix  dT*f  resembles  the 
matrix 

"c  i  O' 

_0  0  I 

and  dT*j.£dT*.}.  resembles 


the  lower  principal  submatrix  of  E.  It  follows  that  for  p(t)  sufficiently 
large,  the  bivariate  random  variable  (CBf-T|(Rf,  Af,  E+)D,  CF+"T^(R-j-,  A^-,  E+)D 
has  as  its  variance-covariance  matrix  a  matrix  resembling 

°22+yII  °23+yI2 

_°32+y2I  °33+y22_ 

where  the  dependence  of  the  matrix  on  the  trajectory  {(p(t),  a(t),  e(t))|teT} 
has  been  nearly  eliminated.  Clearly  the  variance-covariance  matrix 
associated  with  any  other  time  t’  will  also  have  the  above  form  if  p(t’) 

*  p(t).  It  should  be  noted  that  the  estimate  of  r  will  be  dT*+Eddf  *-j--E  ’ , 

where  E’  is  the  lower  principal  2x2  submatrix  of  E.  As  p(t)  increases, 
the  matrix  dT^Edf 1  becomes  increasingly  less  dependent  on  the 
trajectory  and  gives  an  increasingly  more  valid  representation  of  r. 

In  addition,  an  examination  of  dT*+EdT*f  yields  that  terms  of  the  nature 
of 


pflT 

will  be  present,.  Clearly  then  if  it  is  desired  that  these  terms  go  to 
zero,  p(t)  must  be  large  with  respect  to  X.  This  may  be  facilitated 
if  a  configuration  can  be  chosen  such  that  X  is  small  with  respect  to 
ranges  which  may  be  encountered.  If  such  a  configuration  is  not  possible, 
practical  attempts  at  estimating  r  in  this  fashion  will  fail,  and  in 
general  will  result  in  overestimation  of  the  magnitudes  involved.  For 
this  reason,  it  is  clear  that  the  magnitude  of  X  required  for  successful 
evaluation  of  dyncmic  parameters  will  be  so  small  that  curvature  problems 
are  of  no  consequence.  It  should  be  mentioned  at  this  time  that  the 
vectors  of  means  of  the  members  of  the  family  {T*(R+,  A+,  Ef)|teT}  are 
biased  away  from  the  respective  vectors  in  (T*(p(t),  ot(T),  e(t))|teT}. 
These  biases,  however,  may  be  shown  to  be  very  small  when  p(t)  is  large. 
The  demonstration  of  this  fact  is  routine  but  tedious  and  will  be 
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omitted.  It  follows  that  if  p(t)  is  large,  the  bivariate  process 
{(Q8f  -  A^.,  E-j.)],  QFf  -  ^2^t,  ^t»  E+)3)|teT>  will  have  mean 

near  zero  and  independent  of  t  in  T.  In  this  portion  of  the  trajectory 
the  process  is  very  near  to  weak'y  stationary  and  hence  behaves  as  a 
sample.  One  may  then  estimate  the  diagonal  members  of  the  variance- 
covariance  matrix  of  this  process  and  hence  obtain  an  estimate  of  r. 

In  Table  I,  estimates  of  r  are  derived  using  FPS-16  radars  and  AN/GMD-IB 
windsets  at  various  base  lines,  or  X  values.  The  scarcity  of  data  is  re¬ 
cognized,  but  the  tendencies  indicaled  bear  out  the  theoretical  develop¬ 
ment,  and  would  seem  to  indicate  that  the  AN/GMD-IB  is  operating  at  near 
+o  engineering  specifications  insofar  as  measurement  of  azimuth  and 
elevation  angles  is  concerned. 

The  problems  of  greatest  magnitude  associated  with  The  rawinsonde  system 
appear  to  be  related  to  errors  in  height.  A  systematic  study  of  the 
effects  of  height  errors  on  position  estimates  for  the  rawinsonde 
system  is  found,  for  example,  in  DO  and  indicates  That  for  small  eleva¬ 
tion  angles,  system  performance  depends  critically  on  height  errors, 
with  small  excursions  in  height  causing  large  excursions  in  correspond¬ 
ing  position  estimates. 

Errors  in  height  may  be  divided  into  two  classes.  The  first  are  errors 
due  to  systematic  biases  and  random  fluctuations  in  the  instruments  mea¬ 
suring  temperature,  pressure,  and  humidity.  The  second  are  errors  due 
to  physical  situations  which  may  cause  variations  in  the  validity  of 
the  hydrostatic  equation.  It  is  clear  that  errors  such  as  these  can 
vary  with  height  and  from  run  to  run.  If  an  individual  trajectory  is 
examined,  only  the  excursions  due  to  random  fluctuations  in  the  instru¬ 
mentation  of  the  radiosonde  package  are  apparent.  The  remaining  errors 
generally  take  the  form  of  biases,  the  presence  of  which  is  not  obvious 
from  the  examinafion  of  a  single  trajectory  unless  this  trajectory  is 
compared  with  a  trajectory  measured  by  a  device  whose  error  behavior 
results  in  a  trajectory  which  is  not  so  greatly  biased  away  from  the 
absolute;  a  trajectory  such  as  that  measured  by  an  FPS-16.  The  exact 
error  behavior  of  height  measurements  due  to  random  fluctuations  in 
temperature,  pressure,  and  humidity  is  difficult  to  obtain,  even  approx¬ 
imations  of  this  quantity.  This  is  due  to  the  complexity  of  the  mathe¬ 
matical  operations  involved  and  the  number  of  interacting  factors.  How¬ 
ever,  the  portion  of  the  variance  of  height  estimators  due  to  this  aspect 
is  believed  by  the  author  to  be  quite  small  in  relation  to  the  remaining 
contributors.  This  explains  the  smoothness  and  small  excursions  observed 
in  a  single  trajectory.  It  Is  only  when  many  trajectories  are  available 
which  may  be  compared  with  absolute  trajectories  that  the  effects  of 
the  remaining  contributors  become  apparent.  It  is  suspected  by  the 
author  that  the  effects  of  the  latter  contributors  may  be  an  order  of 
magnitude  greater  than  the  former.  As  has  been  remarked,  the  magnitudes 
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RESULTS  OF  DYNAMIC  PARAMETER  ESTIMATION  FOR  TWO  AN/GMD-IB  WINDSETS 
UTILIZING  THREE  FPS-16  RADARS  AT  VARIOUS  SLANT  RANGES. 
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of  the  errors  in  height  may  be  expected  to  vary  with  height,  and  from 
run  to  run;  a  single  error  parameter,  or  for  that  matter  any  number  of 
error  parameters,  may  not  suffice  to  describe  the  system  adequately, 
in  this  case  one  is  forced  either  to  calculate  bounds,  if  possible,  or 
to  seek  situations  for  which  error  parameters  are  stable.  Toward  this 
end  the  following  procedure  may  be  useful  for  determining  average  abso¬ 
lute  height  errors  for  the  raw  insonde  system. 

Let  6(t)  be  defined  for  each  t  in  T  by 

6(t)  =  p(t)sine(t)  +  Xsinn. 

Here  p(t),  e(t),  X,  n  have  the  meanings  ascribed  previously.  It  is 
not  difficult  to  see  that  0(t)  is  the  height  of  the  balloon  under  ob¬ 
servation  at  time  t  in  T  as  seen  from  the  AN/GMD~(  ).  A  reasonable 
estimator  of  this  quantity  based  on  observations  from  the  radar  is 

8^  =  Rj.  sin  E.J.  +  Xsinn. 

For  any  t  in  T,  by  the  same  logic  as  presented  earlier,  the  variance  of 
0.J.  is  approximated  by 

Var(0^.)  %  sin2e(t)Oj|  +  p2(t)cos2e(t)a33. 

Now  suppose  that 


Then 


p2(t) 


or 

p2(t)o33  $  a | |. 

If  It  is  assumed  that  0  <  e(t)  $  n/2,  cosc(t>  >  0  and  p2(t)cos2c(t)o33 
$  cos2e(t)a||j  hence,  p2(t)cos  2e(t)a33  $  ( I— si n2e(t) )a| |  and  p2(1 )cos2e(t)o33 
+  sin2e(t)a||  $  a  1 1 .  From  this  it  follows  that  Var0^.  $  a  j  |  when 
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(Accepted  values  C5D  for  the  error  behavior  of  the  FPS-16  give  a,.  =  15 
feet,  032  "  *01  degree;  hence,  the  erroi  in  8-j-  will  be  bounded  by  i 5 
feet  for  slant  ranges  up  to  83,000  feet.) 


If  Z-j.  represents  the  estimate  of  balloon  height  from  the  rawinsonde 
system,  it  follows  immediately  that  Var(Z^.)  5  Var(Z-}--0j.)  -  a 1 5  as 
long  as 


p(t)  s  . 
033 


Observing  that  one  may  choose  a  portion  of  a  given  trajectory  to  satisfy 
this  condition,  and  also  that  for  any  t  in  the  given  range  the  expected 
value  of  (Z-j--0.j.)  should  be  zero,  an  estimate  of  an  average  variance  of 
(Zf-e-f.)  may  be  calculated  for  this  interval  and  hence  a  lower  bound  for 
Var(Z-j-)  in  that  range*  It  is  again  to  be  stressed  that  such  estimates 
will  vary  from  run  tc  run  and  could  only  be  expected  to  yieid  gross  be- 
havior  characteristics.  However,  in  many  circumstances  gross  behavior 
characteristics  are  better  than  wrong  behavior  characteristics  or  none 
at  all. 


APPENDIX 


A  General  Statistical  Approach 

Lei  X  =  (Xj,  %2>  Xn)  be  a  multivariate  random  variable  with  mean 

x  =  (X| ,  x2,  . ..,  xn)  and  variance-covariance  matrix  E,  Let  T  be  a 
transformation  on  Euclidi  :n  n-space  of  class  C1  in  a  region  R  contain¬ 
ing  x.  Suppose  that  T  has  parametric  representation 


Y |  ~  x2>  **•* 

y2  =  f2(x|  ...  xn) 


y  =  f  (x.,  x„,  ....  x  ) 

•  n  n  1  9  *  n 


Considei  the  multivariate  random  variable  V  =  (/ j,  y2>  *  yn)  obtained 
from  (Xj,  X2,  ...»  X„)  and  T  in  the  obvious  fashion.  From  the  multi¬ 
variate  random  variable  Y  =  T(X)  it  is  desired  to  recover  information 
about  X,  in  particular  the  vector  of  means  and  variance-covariance 
mafrix.  The  following  theorem  will  be  required  which  is  similar  to 
Theorem  25  of  chapter  f  of  C33>  and  is  included  for  reference. 


Theorem  I ;  Let  T  be  a  transformation  on  Euclidian  n-space  of  class  C1 
in  a  region  R  containing  x.  Let  dT  represent  the  differential  [3]  of 
T  at  x,  and  assume  dT  is  nonsingular  or  equivalently,  that  the  Jacobian 
of  T  does  not  vanish  at  x.  If  y  -  T(x),  then  there  exists  a  neighborhood 
J  about  x  which  is  mapped  in  a  t-l  fashion  onSo  a  neighborhood  T(J)  of 
y,  and  on  T(J),  a  transformation  T”  may  be  defined  of  class  C1  in  T(J) 
such  that  for  any  z  in  J,  T"*(T(z))  =  Z,  and  moreover,  if  dT“  is  the 
differential  of  T”'  evaluated  at  y,  dTdT"  =  I,  when  I  is  the  identity 
matrix. 

The  following  theorem  may  now  be  established. 


Theorem  2:  Let  X,  T  be  as  previously  defined.  Let  t  *  be  as  in  theorem 
77  If  TM|  *s  defined  for  <?'  n-vector  v  by 

imi  =  /TTf , 

j=l 

then  as  is  a  measurable  transformation  there  exists  a  number  n>0 
such  that 
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Preceding  page  blank 


Pr{X^T~ 1  (JO  }  5  • 

n 

Here  Pr  denotes  probability,  tr(E)  the  trace  of  E. 


Proof:  Choose  n  so  small  that  {z||  z-x||<n)  is  contained  in  J.  This 
assures  that  T“  is  delined  on  ?({z  | | z— x | | <n } ) ,  hence  {z| | lz~xj |<n}  is 
contained  in  {z|z=T“* (T(z))}.  Now  consider  the  event  {X^T”! vT)}. 

(Observe  that  T“*  will  be  measurable.)  It  is  clear  that  {Xj*T~  (lO}e{  |  |X-x|  |^n) 
so  that  PKX/T”  (W}  $  Pr{  |  |X-x|  l^n).  However, 

n 

Pr{ | |X-x| l^n)  =  r{  E  (X-x)2  5  n2}. 
j=l 

From  this  it  follows  (see  [6],  section  15.7)  that 


n 


Pr{ | |X— x | |^n) 


$ 


Ea . . 

j=lJJ  _  tr(E) 
^ - 


This  proves  the  theorem. 

This  theorem  states  that  the  probability  that  X  and  T  *00  will  differ 
depends  on  the  size  of  the  trace  of  E,  and  on  the  area  of  definition 
of  T"  .  If  the  nature  of  T  is  such  that  n  is  very  small,  unless  tr(E) 
is  correspondingly  small,  it  may  occur  that  in  a  probability  sense  V  may 
furnish  little  information  about  X. 


Estimation  of  Parameters 


To  estimate  the  matrix  E  from  a  sample  of  the  multivariate  random  variable 
j,  one  may  choose  to  estimate  the  variance -covariance  matrix  of  V,  call 
it  A,  and  then  take  dT^Adf-*  as  tne  desired  estimate.  This  estimate 
will  be  biased,  especially  since  the  expected  value  of  y  is  not  necessarily 
T(X).  However,  if  tr(E)  is  sufficiently  small,  this  will  not  be  signi¬ 
ficant.  In  dealing  with  multivariate  stochastic  processes,  affairs  be¬ 
come  more  complicated.  Let  {X-j-|teT}  be  such  a  process.  Let  the  vector 
of  means  of  X^.  be  givan  for  any  time  t  in  T  by  the  vector  function  X(t), 
and  let  E  be  the  variance-covariance  matrix,  independent  of  t.  Let  T  be 
a  measurable  transformation  on  Euclidian  n-space  of  class  C1  in  a  region 
R  containing  {X(t)jteT}.  For  each  t  in  T,  let  y +  =  T(Xf).  The  multi¬ 
variate  process  {^flteT}  so  defined  will  no  longer  have  variance-covariance 
matrices  independent  of  T.  Let  A-f-  denote  the  variance-covariance  matrix 
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I 


I 


of  Y-j-  for  each  t  in  T,  Observe  that  no  nontrivial  estimate  of  Af  is 
possible  from  a  single  realization  of  {^|teT},  hence  a  technique  such 
as  described  above  for  finding  Z  is  not  possible. 

$ 

It  has  been  shown  C2],  however,  that  for  each  t  in  T  there  exists  a  trans¬ 
formation  Lf  on  n-space  such  that  Lf(Xf)  is  a  multivariate  random  vari¬ 
able  with  mean  T(X(t))  and  variance-invariance  matrix  dTfZdTf,  where 
dTf  represents  the  differential  of  1  at  x(t).  Moreover,  for  any  chosen 
precision  figure  e,  there  exists  a  number  6(e)  i  uch  lhat 

i 

Pr{||y+-  Lt(Xt)|U  e>  . 

! 

If  for  a  chosen  e>0,  tr(Z)  is  sufficiently  sma:l  I  that  (L-j-(Xf)  |teT}  is 
a  uniformly  satisfactory  substitute  for  {y+|teT>,  the  following  approach 
is  possible.  For  i  and  j  among  I,  2,  ...,  n,Jet  y(i,j,t)  and  a’(i„j,t) 
represent  the  i,  jth  entries  of  dTf  and  dTfZdTf,  respectively.  Nov/  . 
assume  Z  is  diagonal  so  that  one  may  always  write 

n  i  i  • 

o'(i,j,t)  =  Z  y(i,k,t)y(j,k,t)o. .  .  (I) 

k=l  kk  (  , 

Since  dTfZdTf  is  symmetric,  it  is  possible  to  find, 
n(n+l ) 


distinct  equations  of  the  above  form  of  which  at  most  n  may  be  independent. 
A  reasonable  choice  of  n  equations  is  obtained  by  setting  i  =  j,  obtaining 
for  Z  =  I ,  2,  . . .,  n. 


a'(Z,Z)  -  Zy2U,k,t)a, .  . 

i _ i  kk 


Observe  that  the  fact  that  dTf  is  nonsingular  does  not  Imply  that  this  > 
system  is  nonsingular.  Now  suppose  that  for  Ct|  an  interval  in  T, 
a(i,j,t)  and  y(i,j,t)  are  integrable  functions  of  t.  One  may  then  write  t 

from  ( I )  i 


x~r-  J  o’(  i  ,j,t)dt 
T2"TI  Tl 


Z  C+-4-  /y(i,k,t)y(j,k,t)dt]a 

k--|  VT I  t, 


I 


kk' 


(2) 


J 
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tut  ifmmta* 


T 


IN»f* 


I 


I  *  ; 

•Like  (I),  (2)  is  also  a  linear  equation  'in  the  okfc  with  the  coefficients 
completely  determined  by  the  knowledge  of  {x(t) |teT}‘ and  t.  If  the 
integral  :  1  ! 


f  '  <  1 

+  4+-  /2ff'Ci.j»+)dt  '• 

*  1  T2  T!  ft,  •  •  :  »  *  ^  f 

is  known,  for  some  n  distinct  values!of  (i,j).,  and  the  resulting  system 
is  nonsingular,  the  arfe  completely  determined.  However^  these 
values  are  not  known.  If  es+imates  of  some  satisfactory  set  of  these 
integrals  exist, i however,  this  may  lead  tq  estimates  of  the  o^k. 

i 

Consider  the  system  defined,  far  £=1,2,  ...,;n  by  t 


i 


t  1  ‘ 

_J _  ,2  n  i  |  +o 

t,-t.  ;  b*(j,j,t)dt  =  z  /V(j,j,t)dt:a 

'*  Tl  .  k=l  T2  Tl  t 

4  I  .  , 


Let  [t|,f2]  be(  divided  into  m-l  equal  intervals  of  length  A  by 

t|  =*S|  <  s2  <  ...  <  sm  =  t2.  Let  X-  =  Xsj  for;j=l,  2,  ...,  m-and  let 

Xj  be  the  mean  of  X,..  Define  {L  j(Xj')  |  j=s|  ,m}  in  the  obvious  manner,  and 

let  L:  j  be  the  ithJ component  of'Lj.  Now  for  <i  =  I,  2,  ...,  n,  define 
J,i  ,  J  i 

•  i  .  Vn 


'  i 


m- 


(X.) 
J  . 


T.Uj)]2  . 


i 


It  !fol lows  that 


.  m  :  .  m 

?  Ei(V  "  W32  » 


.  ’  m 

■  ■ v  •  Z  o’ ( i , i ,s . ) A. 
* I  •_!  J 


2  *1  j=l 

J  • 

Since1  a’(i,i,t)  is  integrable  on  [tj,t23,  as  A  approaches  zero, 


m 


-~r~  2  o»(  i,i,s.)A 

t2“tI  j=l  J 


? 


1 
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i 


i 


approaches 


V+| 


i » 


t)dt. 


The  I  i  near  system 

n  .  t9 

S.  =  E  Cq—r  /  W2< j,k,t)dt3o.  j=l,  2,  n 

J  k=l  +2~+l  t,  k 

may  now  be  solved  for  providing  that  the  system  is  nonsingular. 

If  it  should  occur  that  the  system  is  singular,  other  members  of  the 

n(m+l ) 

2 

distinct  equations  of  this  form  may  be  tried  or  Lt|,t2D  varied  to  obtain 
a  sufficient  number  of  independent  equations.  ! 
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